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ABSTRACT 


This investigation analyzes the operating character- 
istics of an ejector having a heated primary jet and a 
constant area mixing section. The effects of several key 
parameters are determined on the basis of idealized theory. 
Performance characteristics are described in terms of certain 
non-dimensional coefficients obtained by a systematic dimen- 
Sional analysis. 

In addition, experimental results obtained from a model 
ejector are presented. These results are in reasonable 
agreement with the theoretical analysis. 

Suggestions are made for iuruher improvements bork 


in the theory and in the experimental facility. 
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ERRATA 


Dade uepsetimen 7 Should read "tieorerical’™ vice 
Mimeoreteial " 


mage 25 pata 2.2.a, line 2yshoullastead “suitably” 
Vice sumtable!! 


made 15, para 222,50, last Pune should read “suitably” 
vice "suitable" 


page 22, line 7 should read"downstream" vice "down stream" 


page 28, lines 1 and 2 should hyphenate the word 
"performance" between the m and the a 


page 62, third line from bottom should read "of an" 
Vices on an” 


age 68, eighth line from bottom should read "iteration" 
pag 
VEeCeomuiativeratlon 


page 87, item Tr,/TTy is described as 


"Average exit total temperatures.", this statement 
Should read "Ratio of the secondary-to-primary total 
temperatures," 


List of Symbols and throughout the thesis where Cp 


has been used as the symbol for specific heat 
at camstant pressure, Cp should be changed to read 


“ye 
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Pie |eCcuoCt., or Jel DUMP, 1S 4 GCevice wate uses the 
relatively high energy of a small heated jet to produce a 
pressure rise in a relatively large mass flow. Such devices 
have many diversified applications. Jet pumps can be used 
as boundary layer control systems to delay the effects of 
ibew Separations. Air may be either plown over the lifting 
surface or sucked away from it. Similar devices, greatly 
enlarged, have also been employed as jet engine test cells 
where a large mass flow at high velocities is required. 
(Ref. 3). Further applications are in the fields of thrust 
augmentation and STOL/VTOL aircraft. (Ref. 6 and 7). 

It was desired to study the effects of several key 
parameters upon which the performance of such ejectors depend. 
It became evident that before a specific parametric analysis 
could be performed, some meaningful set of performance coef- 
ficients should be developed which describe adequately the 
operation of an ejector. Such coefficients would also aid 
mm cbhe design of j,e%6 pumps in @enerat er ics soewioriance 
coefficients have been derived on the basis of a logical 
dimensional analysis of the principal parameters which control 
the operation of such devices. 

Once the coefficients were defined, a systematic study 
was performed to determine the effects which several key 


parameters have upon ejector operation. 


This study is based upon the continuity, energy, and 
momentum equations as applied to the axi-symmetric, two- 
dimensional mixing of two concentric jets in a constant area 
duct. The final numerical results of this parametric angie au. 
were determined by means of the Control Data Corporation 
1604 digital computer. It should be noted that losses 
associated with friction and with incomplete mixing have 
been neglected, so that the results obtained represent the 
optimum performance levels which an ejector may approach but 
never exceed. 

An experimental ejector, similar to that presented in 
Reference 2, was fabricated and installed in Building 216 at 
the Naval Postgraduate School, Monterey, California. A 
brief test program was conducted with this device in order 
to obtain an experimental comparison with the idealized 
theoretical trends. 

Both the theoretcial analysis and the experimental test 
program are extensions of work anitdaved by Belvermein 
Reference 1 and continued by Ridder and Summers in Reterenee 
2. In the case of the parametric analysis, however, some 
results obtained in Reference 2 were found to be in error. 
Specifically, the tabulated mass flow ratios shown in 
Appendix C of this reference do not agree with values obtained 
by subsequent calculations. For this reason, along with the 
desire to develop a more simplified approach, the theoretical 
procedure presented in this report is basically original. 

The results obtained in Reference 2 apply only to compressible 


flows, while the methods presented herein may be applied to 


da@ 


either compressible or incompressible flows. The specific 
design of the experimental ejector suggested in Reference 2 
was scaled downgto approximately one-fourth the original 
linear size to meet cost and power supply limitations. Also, 
the actual ejector tested is of the variable inlet pressure 
type as compared to the variable back pressure system 
originally proposed. 

Thereteore, the goal of this Studvenas. Lo dev clon soamere 
basic set of performance coefficients which adequately 
describe the operation of an ejector; to investigate the 
efirect which several key parameters @aeve Upon these perfor— 
mance coefficients; and to obtain a comparison between 
experimental results and the idealised Cheomm 

The author wishes to express his sincere appreciation 
for the invaluable assistance and guidance provided by 
Professor Theodore H. Gawain ofethe Naval Postgraduate 


school, Monterey, California. 


2. THEORETICAL ANALYSIS 
2.1 Dimensional Analysis 
It may be seen from Figure 1 that the main parameters 
governing this problem are those which define the thermo- 


il 


@ynamie states of the primary fluid P Ti )3 the 


ale: Pi Tale 
secondary fluid (Pings 72 Tino » T,)3 and the fluid at the 
exit from the ejector (Pin » Pos Tin » T2)- The amount of 
mass in each system and the geometric relation between 


the individual systems are also of great importance. For a 


Chermodynamic system. in which electriveal Citeers ave 
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insignificant, four basic dimensions suirfice Co nen—-dimen— 
sionalize all possible physical characteristics of tne system 
For this study, the basic dimensional system of force, 
length, time, and temperature has been chosen. It is felt 
that by judicious selection of four reference parameters, 
which collectively represent the four basic dimensions, a 
more meaningful set of non-dimensional performance coef- 
ficients can be defined. 

After analyzing the overall operation of an ejector, 
the following four parameters were chosen to represent the 


basic dimensional system: 


Symbol Unit@s 
1. The ambient density of the cm PTS 
secondary fluid 1A 
2. The ambient temperature Uno 0 
of the secondary fluid 
3eeelne mass flow raver orm my Pay 
primary fluid L 
4, The available enthalpy drop H., Le 
2 
a 


Using these as a dimensional base, it is possible To 
express any other physical quantity in the form of a non= 
dimensionalized coefficient. (See Buckingham-Pi Theorem, 
Appendix A). 

c.e Definition of Dimensionles eer ciCnys 
a. Pumping Power Coefficient 
The pumping power coefficient is defined as the 
maximum available power at the exit from the constanveaee 
mixing section, suitable non-dimensionalized by means of the 


dimensional reference parameters mentioned above. 


im 


._ (m, + m,)H 
Cn il ae ir) 
H 


my a 
b. Suction Power Coefficient 

An ejector is capable of inducing a flow of secondary 
fluid from some pressure below that of static pressure at 
the exit from the mixing section. From a thermodynamic 
analysis (Appendix A) it can be shown that there is some 
ideal equivalent power required for this type of isothermal 
compression. The suction power coefficient is then defined 


as this equivalent power, suitable non-dimensionalized. 


co 06e Mo Tyo AS 


7" ene (2) 


m Ha 
Ge. tloval Power Coefficient 
With an ejector developing power both in the form 
of available exit enthalpy and in the form of isothermal 
compression of the secondary fluid, it is appropriate To 
define the total power coefficient merely as the sum of 
these two outputs. The coefficient so established represents 


the maximum possible power developed by an ejector. Thus 


= + 
eae ome Ds (3) 


ad. Momentum Coefficient 
The momentum coefficient is defined as the output 
momentum at the exit from the constant area mixing section, 


properly non-dimensionalized. 


3 


It should be noted that the power coefficients Cop and COs: 
being the ratios of output blowing and suction power, 
respectively, to available power, really represent the 
efficiencies of the system as a pump and as a suction device. 
The total power coefficient therefore represents the over- 
all efiwedsency : 
See Driging PressurésCoeiiikedicrad 
The difference between the primary supply pressure 
(Pin? and the static pressure at the exit plane (P.) is 
non-dimensionalized by a suitable reference pressure. This 
then defines the driving pressure coefficient. 
oe fh (5) 
Ppa Ha 
f. suction Pressure Coetiiement 
A second pressure parameter must also be defined 
for an ejector with a variable inlet (suction) pressure. As 
before, the difference between the secondary fluid total 
pressure (Pro) and the exit static pressure (P,) is properly 


non-dimensionalized to define the suction pressure coefficient. 


The major portion of Ghrs resort deals Waves 


compressible operating fluid and it became more converiiems 


P2 Pmo : 
tO use pressure ratios, and ==, as controlling pa 
Pry Pry 
meters. There is, of course, a correlation between tia. 


ratios and the above pressure parameters as shown in Figure 
28. Hence, the results of the study can be applied to either 


the compressible or incompressible situations. 
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g. Mass Flow Ratio 
The rave of mass flow induced through the secondary 
Syovem Dy the ac@on of the primary jeu is of “primary 
iiporvance 1n Gescriming the pumping Gharacveristics of an 
ejector. A non-dimensional mass flow is therefore defined 
as the secondary mass flow rate (mo) divided by the primary 


or supply mass flow rate. 


X = m2 Cp 


h. Area Ratio 
Trre geometrical factor which controls the performance 
of a jet pump is the physical sizé of the conponent systems. 
mere is, therefore, a need to define certain Size parameters 
Mmenmconmplete the non-dimensional des¢rapumoeneor al eleevom 
syovem. AS Shown in Appendix Al Gneverecs .secrtonalewawea Of 
Enew primary Jel mozzle, and’ the COmstamesarca Mixmme seevion 


ea be nNon—dimenstonalized to give twWomareameeehutetcnps. 








oie 5 | Ha (8) 
eal 
and 
a Pmoyy 
eh. A a ae (9) 
ay 


It is evident, however, that the individual size of 
each component is relatively unimportant and that the 
controlling geometric factor is simply the ratio of the 
mixing section area (Az) to the primary nozzle area (A,). 
For this reason, the above defined area coefficients were 


replaced by the ratio A2/A, and Ca: 


i? 


With these coefficients, it @e@eecciblesto decerine 
thewessential operating characteristagsmerean C)}ecver syverenm- 
This description is essentially independent wot the. operating 
fluid or of any absolute dimensional magnitudes. This method 
should aid in the specification of ejector design requirenecnrs 
2.3 Theoretical Development 

This analysis is an investigation of the physical 
relationships which govern thesoperatvioneof a ject \eunipe 
The resulting equations may be solved by a systematic itera- 
tion technique which can be performed on any digital computer. 

Figure 1 shows the station locations and symbolic 
notation employed throughout this development. 

The three basic conditions which must be satisficedigs 
the constant area mixing of two streams of fluid in an axi- 
symmetric duct are the equations of continuity, momentum, 
and energy. These equations, applied to the frictiomlea. 


flow with uniform velocity distri ons, are. 


Continuity m4 + m. - m (10) 
Momentum PoAy + PoAo — Poh, = n3V. = m,V5~m,V, Cine 
Energy miCpTpz + m3C,Tao = m3C,Tp2 Gi 


The detailed solution to these equations iS given wm 
Appendix A, and will not be covered here. The following is 
a summary of the resulting equations obtained, along with a 
listing of pertinent parameters. It is felt that in this 


manner the method of solution can more readily be shown: 
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Basic Equations 








P y-1 =i 
Peyieri \ y ews 4 1/2 ote 
i ‘ toe M5 i | 
. Yel 1/y-1 
ae Tp 3 1+ om 
aon — t.. |) | =e (14) 
: T1/\*T) me M2 
th M5 \ fir ye" wn 21172 
Keene) {| ° |/| 222) ee (15) 
| \\ tre || lee 
my : 1+ "3 Mp 
P2-Po 
hs Pov, 
Or 
d,s {_i_ P3 ee a 1+ Gm? 
4 SyM.-} \{P 
Vogal Tl an , patele 2 2 
ig 
Gilg) 
2 
9 = 1 Sl, a a 
M x2 LS opv.2 32(Ke-1) ae. 
a = 6 K 
where 
a 
AP, ; + Y-1y : P ‘ 
5 = 2 528 a v¥/Y~- 
e 2 
oer al oy, , + Y= My 
iE i 
SI een ee Ta Gees) 
a= =| = + X 
Tipo Tipo 
ie 
P 
—- +O) ym 
23 | PoVq x2 1 (19) 
om ios iE a a 
; — : ia 1 mo° 


—~ bie 
Lt t= ———_—_—_—— -— X (20) 
a ) 
me /.0 


Known Parameters 


> 


Method of (selurion 


Assume a trial value for Ms: 


Assume an initial trial value ™or pf equal to unas 


Calculate M 04 kes D es and © 


qe 
Define an error function 


M" 


ERF = 9, - © 


M so 


Vary the trial value of M. until ERF = 0 for initial 


2 
value of u. 

Using equations 6 thru 8 calculated a new value for u. 
With the new u start Ms iteration again. 

Repeat until two successive values of M5 are equal. 

This gives the proper value of M5 and u for the specified 
constant parameters. 

Once M. and uw have been obtained, the performance 
coefficients may be calculated. 
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These computations are performed by the Fortran 


program JETPUMPI as shown in Appendix A. 


See APE RI MENIAL TST Pacman ry 

The model ejector employed for this study is essentially 
a scaled-down version of the design originally proposed by 
Ridder and Summers in Reference 2. A scale factor of 
approximately one-fourth was necessary to meet financial as 
well as power supply limitations; further, minor modifications 
were also made to the paisi cl deeteem 

The original proposal in Reference 2 called for a 
control valve to be installed at the exit from the mixing 
BecuLon.  suechea valve would allow foumcontrollinge the 
static pressure at the exit to the ejector. Such a device 
is called a variable back pressure ejector. The actual test 
ibacllity 1S Open vo avumospheric Pressmieoeameenec Cxay plate 
and therefore not of the variable back pressure type. A 
control valve is installed in the secondary air system to 
facilitace Convrolling the Secondary sameeren me DCe sours. 
The pressure can be lowered from atmospheric to a partial 
vacuum. For this reason, the ejector tested is of the 
Variable suction Uype. 

The test model is comprised of three basic systems: 
The primary or high energy air system; the secondary or low 
energy system; and a variable length constant area mixing 
Seco iar 
3.1 Primary Air System 

The primary air system, shown in Figure 2501s the 


energy source for the ejector. High pressure air is supplied 


Ie 


from a Pennsylvania Pump and Compressor Company HAE compres- 
sor and is delivered via a convergent nozzle Vo the mixing 
section. The total pressure, or pressure energy, supplied 

to the primary system can be controlled by a Cash-Acme pressure 
reducing and regulating valve assembly just upstream of the 
primary plenum chamber. In the present tests, the ejector 

was operated over a range of 4 to 25 psig although the 
compressor is capable of supplying air at pressures up to 

OOM ersare. 

The plenum chamber acts as a settling tank to minimize 
any turbulence or fluctuations in the primary flow. The 
air passes from the plenum through a bell-mouth into a two- 
inch pipe. Installed in the pipe is a stainless steel 
square edged orifice plate along with four flange static 
pressure taps and four vena contracta static pressure taps. 
An iron-constantan total temperature probe is also installed 
just upstream of the orifice plate. This arrangement makes 
it possible to measure flow rates in accordance with ASME 
standards. (Ref. 5). 

Downstream of the orifice is a continuous burning com- 
bustion section. (Figure 3). This combustion section is 
fabricated of 304 stainless steel two-inch tubing. The 
actual combustion chamber is located in the middie of “Gils 
section and a butterfly valve has been incorporated to con- 
trol the amount of primary air passing through the combustion 
chamber. The remaining primary air by-passes the chamber, 
and in this manner the total temperature, or heat energy, 


of the primary @itrecan be copmtcorlica- 
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After leaving the combustion section, the high energy 
air discharges into the stainless steel mixing section 
through a convergent nozzle. The nozzle, as depicted in 
Figure 4, is located inside the secondary air system plenum 
chamber in such a way as to develop two concentric jets at 
ene entrance Lo tne mixing sect item: 

3.2 Secondary Air System 

The high energy jet issuing from the primary nozzle 
imeees a tlow of Secondary alr, fromeamprenuy condivLoma, 
through a bell mouth inlet into a three-inch standard pipe. 
Mm square edged orifice system is installed in the pipe for 
mass flow rate determination. The system is identical to 
that employed for the primary air. 

Downstream of the mass flow rate measuring station, 
emere is a three=inch sate valve. Syewarying tne seouum 
or Clic Valve, une total pressure Of sume oe comaary cil mmecan 
be decreased below atmospheric. In this manner the opera- 
Dion Ol va SUCGlOm type €jJ/eCClCOr cans bemcamulraved. 

[The secondary air then envers a ptenumechanber, vne 
purpose of which is to provide more nearly uniform and axi- 
symmetric inlet conditions. The air then enters the mixing 
section through a convergent nozzle. 

S25 Mixing secon 

Both the primary and secondary jets meet at the entrance 
GO a Constant area duct or mixing Section. Mixing secures. 
of three different diameters are available. Also, primary 


nozZles of three different sizes can be installed. Therefore, 


ea 


nine different area ratios can be Obtained =) 'or Gotomeuua, = 
however, only one mixing section was used. 

The section employed is made of nominal two-inch 
stainless steel tubing (ID = 2.25 inches), fabricated in 
three sections. By varying the combination of sections, 

a range of L/D ratios from two to ten can be obtained. L is 
the length of tubing down stream of the nozzle exit. 

Five instrumented stations (Figure 5) are equally 
spaced along the length of mixing section. At each station 
are two averaging static pressure wall taps (1/32" diameter). 
There is also a 1/8-inch hole with a boss assembly for 
obtaining pitot static and total temperature surveys at 
Gach Station: 

For normal applications, a diffuser would be installed 
downstream of the mixing section to convert any velocity at 
the exit to a static pressure rise. For this study, NOwWeraae 
a diffuser was not incorporated since the mixing process 
was of prime interest. 

3.4 Instrumentation 

All pressures were obtained from a 94-inch mercury 
manometer bank with the exception of the AP across each 
orifice plate and the exit traverse pressures. The AP's 
were read from a water micromanometer calibrated in Cemeame 
meters, while the traverse pressure was read in inches of 
water from a 48-inch single tube manometer. The two tempera- 
tures upstream of the orifice plates and the exit total 
temperature were measured with iron-constantan thermocouples 


with a Leeds and Northrup millivolt readout. 
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The mass flow rates were calculated based upon Reference 
4, as shown in Appendix B. The total temperatures recorded 
upstream of the orifice plates were assumed to be the 
static temperature for these calculations. Although the 
mass flow measuring system is equipped with both flange and 
vena contracta static pressure taps, the pressures measured 
with the flange taps were used in the determination of the 
mass flow rates. 

The overall test ejector and instrumentation system 
are shown on Figures 6, 7, and 8. 

Be oLest. Program 

Once the model ejector was fabricated and installed 
in Bldg. 216 at the Naval Postgraduate School, the system 
iMeGerwenuy Severaletrial runs in order to eliminate all air 
leaks and to check out the instrumentation system. Upon 
completion of these trial runs, nine runs were made to 
collect the necessary raw data. 

Mige first =toureruns were made wien eaheag ra ti© sol 
10.125:1 while the last five were with an area ratio of 
20.25:1. All runs were with a temperature ratio of unity 
since the combustion section was not in operation due to 
The Wack ol fueleand a Soni cee iMelyit low orifice: 

A typical run consisted of the following steps: 

1. Establish a given P2/Pmy Tenors 

2. Set the secondary control valve in the full open 

POS@U ton: 
3. Record all data, including an eigh> Sratwon, tea erse 


for both pitot-static and total temperature profiles. 


eS 


4, Close the secondary control valve in four steps 
and record the new data at each step. 
Forsan areasratio of 10.125: IM wheumediffierent vaies 
of P3/Pmy were used. These were values of 0.4, 0.5, 0.6, 
ando0a8. Form the 20.25:1 rationvaltiesmone0.3,0.45055e 
and 0.6 were used. In addition, two additional runs were 
made so that a check could be made upon the repeatability 


of the system. These runs were at an A2/A, of 10.125 and 


5 
P3/Ppy of 0.4.and 0.5. The resulits)of these renuns age 
shewn in Figures 16, 17, L@s@anceer 

The raw data was then reduced as outlined in Appendix 
B. The final numerical results were obtained by means of 


the IBM 360 digital computer through the Fortran program 


EJECTOR. 


4, DISCUSSEON Onan wes 

Figures 9 thru ll show the ideal performance limits 
for an ejector with zero suction, all output being iio 
form of blowing. It was found that a slight decrease in 
secondary total pressure below ambient resulted in a lame 
decrease in mass flow ratio. Also J av zero suction, the 
total power coefficient, and hencejytne overall etficieme 
were Maximum. lor this reason@eemece CUrVes Ceplc pasa 
maximum possible performance levels for this study. As 
previously mentioned, these figures represent the optimum 
performance levels which an ejector may approach but 


never exceed. 
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Figure 9 is an overall mapping of the pumping power 
performance. It should be pointed out that although the 
value of oo decreases with increases in temperature ratio 
(Tiny/Tmo) » the value of the absolute power will increase. 

The apparent indicated decrease is actually a loss in 
Sriere Lemey . 

SIN ia eure LO Sumner Zeomencere lar onship 
between mass flow ratio and the pressure coefficient. It 
ieeseen that C tends to decrease more or less sharply as the 
tiesto {low ratio vs increased. Thigeis in contrast to the 
mewer wrends where the effect of ine€reasingmmass Flow catia 
if nol qusrve so severe. It 18S G@lSercr wmtcrese to Nove that 
(Chnaaorven veMmperavure ratio, the wade wo. Cy is pracrrcadaly 
independent of area ratio. 

The limiting values for output Monenvimn coer! rerent are 
plotted in Figure 11, with mass flow ratio as the independent 
variable. This figure shows that an increase in temperature 
ratio has little effect upon the output momentum, Ca An 
increase in area ratio will, however, result in an increase 
in C,,. 

The detailed comparison between the theoretical trends 
and the experimental results is given im Ficures 2 thru 

ef. The graphs depict the performance of an ejector, boun as 
a power producing device (C and © 


pp ps?» 
and C.). Each set of curves is for a given value of P3/Piny 


and as a pump (Cy 


and A3/A, . The parameter C is considered to be of secondary 


A3 
importance, and detailed values are therefore icrs re. 17-4, 


ao 


however, the range of C values was from 225 to 475. All 


A3 
curves are for a temperature ratio of Minit, since theseen= 
bustion chamber was not in operation for these tests. 

It is readily apparent that the experimental results 
follow the theoretical trends although, as expected, there 
is considerable deviation in actual magnitudes. The power 
curves or and om” show that the greatest deviation from 
the ideal case is a result of a decrease in attainable flow 
rates. This can be attributed to friction and mixinggiie2. 
in the experimental system and also to inlet losses. These 
effects are shown as AX friction and AX inlet, respectijeaee 
The inlet losses result from a decrease in secondary total 
pressure. This drop in total pressure makes it impossible 
to obtain true zero suction operation, and hence, decreases 
the performance level. 

The pumping characteristic curves also show the effect 
of friction, mixing, and inlet losses. The approximately 
constant deviation between the actual and theoretical values 
of suction pressure coefficient is the result of fricraied 
and mixing effects. It reduces the maximum flow capability 
by the amount AX friction as indicated on the diagrams. The 
remaining loss in mass flow ratio is associated with flow 
resistance thru the inlet valve of the secondary system, 
and is denoted by AX inlet. Since the controllable parameter 
for each test was the pressure ratio P2/Pmy s the agreement 
between the theoretical and experimental values of C, is 


D 


necessarily ceiece. 


Dia 


i seneral it can be said that the analytical approach 
employed for this study yields a reasonable overall picture 
of ejector performance trends. The experimental results are 
in agreement with the parametric analysis in regard to the 
qualitative effect of several key parameters upon ejector 
performance. The agreement between theory and experimental 
results is encouraging; however, a somewhat more detailed 
analysis which includes friction and mixing effects will 
mecatily be necessary before actual perrormance levels can 


mecurabely ihe piredicled. 


5. RECOMMENDED MODIFICATIONS TOM UES AC Th i hy 

The experimental model has proven to be useful and 
Pere ble Tor studying ejector perfommames. lt as felt 
uhat the losses shown for this model can be reduced with 
wne LTOllowing modifications: 

a. Install a more efficient valve in the secondary 
air system to reduce the pressure loss through the valve. 

b. Increase the size of the secondary plenum chamber, 
or add some form of baffle system to provide a more uniform 
velocity distribution at the entrance To Che mixing Section. 
This requirement became apparent at high values of Co where 
the flow at the exit plane was not axi-symmetric. 

For an actual ejector, the inlet losses will be reduced 
since the length of plumbing associated with the flow rate 


measuring system will not be required. 


an 


6. RECOMMENDATIONS FOR FURTHER STUDY 

“In order to develop a method whereby the actual perfor- 
mance levels of an ejector system can be adequately 
predicted, the following extensions to present work are 
strongly recommended: 

a. Extend the present theory to include the effecgs 
of friction and mixing losses. This can be accomplished 
by the addition of a friction term to the momentum equation 
and by developing a suitable shape factor to better describe 
TAemMmex tiie pir OCe ssn 

b. Complete the incorporation of a fuel system in 
the existing ejector. By installing a fuel flow regulating 
valve and flow meter, the existing facility will be capable 
of producing a wide range of temperature ratios. 

ec. Conduct a more detailed experimental test progran 
in which velocity traverses are taken at the five Staviem. 
along the mixing section. Such agsvuudy will providers 
better understanding of the mixing process and shoulda yaa 
in the determination of the shape factor previously mentioned: 

qd. Investigate the actual application of ejector 
systems to boundary layer control work. Such an investieae 
tion could be made utilizing the ex@eving 2) ecreremm 
conjunction with a series of aerofolls. By installing ie 
aerofoils in a wind tunnel the effect of the ejector sucvion 


or blowing upon Tlow separattommeould be slpucicde 
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APPENDIX A 


Theoretical Analysis 


1. Dimensional Analysis 

For a thermodynamic system, in which electrical effects 
are insignificant, there are fourwbasre dimensions which 
describe the system. According to the Buckingham Pi Theoren, 
these four dimensions serve as a Dase by which all physieum 
characteristics can be non-dimensionalized. 

Cheosing force, length, timeweand tempecraturesas. cm. 
base dimensions, the following physical parameters were 


ehosen to collectively répresentethese base dimenstonce 


a. secondary ambient density -- pmo 
b. Secondary ambient temperature -- Tmo 
ec. Primary mass flow rate -- mj 


d. Primary available enthalpy -- oe 

With these reference parameters, it is possi plewne 
non-dimensionalize any physical quantity in a thermodynamic 
Sysvem. 

Consider, for a given output power P, a non-dimensional 


DOWe ie © oleic he pC obtained in the following manner, 


oa 
Power - unite som ae 


Therefore, to obtain a power coefficient basca upen 


the selecvedsreterence paraner a]. 


_ PTC 
Pmo a 
if 
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2 
H - Le 
a To 
: FT 
"1 ; a 
Define 
| (H 
C 2 I 
p a bo rc! (A=—1) 


Pp2 Tipo Hg my 
poceings Up fOur equations to sOlveomr emma, 0, c¢, a. 


C = 








FL/T 
. 2\a : : C d a! 
Es 0 LT Ey 
os me i? 
L 
F: atde=l (A-3) 
ice -Ha + 2 -del (A-4) 
oe 2a - 2, +d = -1 (A-5) 
QO: pe 0 (A-6) 
adding equations A-4 and A-5 
-2, =- 0 
or 
(A-7) 
a= 0 
therefore, from equation A-3 
dae (A-8) 
now from equation A-5 
-2a +1:= -l 
or 
Cc = 1 (A-9) 
cee OWET 
ep 
m, He 


in this manner it 1S possible ze define any number wer 


dimensionless coefficients. 


ce. Definition of Performance Coefficients 


Referring to Figure 1, the following parameters are 


obtained: 
= a —_ 
ne = Co (Tiny qT, ) ( A=ieey 
le an 
4 = 30 lem 
H, = C, Tp, (1 = , (A-12) 
oe 
Sa eave 
P2 Voes 
H, = C Tr i {= y (A-13) 
Ping 


From the lst and end Laws of Thermodynamics, 


dh = TaS + V (A-14) 
| dp 
for the isothermal compression from state « to state 
Ty 
an aaa (A-15) 


Ee irei@ ier 


ds . =-v 
4>T2 7 dpy.T2 (A-16) 


for a perfect gas, the incremental chance in entropy mia 
AO) Teweammee Cefined ag. 
AS 


P a P 

4oT2 = =") "T2 = Ry “3 (A-17) 
— po 
Pg NZ 
te 

The equivalent work done by the isothermal compression 


is therefore 


_ P 
a = RTmo gn = (A-18) 
iz 
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Withee se enUnalpics demumedeeic ds now possible to 
develop the power output coefficients for an ejector. 

The power available at the exit from the mixing section 
is defined as the pumping power and therefore the pumping 


power coefficient can be defined as: 


co = (hy + i) (4-19) 
my He 


or substituting equations A-1l2 and A-13, 
mo (E43) (Ss) P-ty’tal 
i ee (A-20) 
1 = [P3/Pp1| a 
The equivalent work done during the isothermal compres- 


sion of the secondary fluid is defined as the suction power, 


and therefore the suction power coefficient is defined as: 


C = Mo “sz, (A-21) 





from equations A-12 and A-18 


Cos = 


— 


zt 
my 


Yr+ 
Tn = (P2/P ny) 5 





| [wre £n P3/Pmno (hea) 








Leesa Uae 2 is defined as the non-dimensional mass flow 
a) 
ratio, x. It is also possible to suitably non-dimensionalize 


the output momentum to obtain the output momentum coefficitenc. 


Ong = il 3 (A-23) 


Or 
fs 
e. (+ xy 23 
H. 
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v= Lene 
aes liga 9) 3] i a aa (A-24) 
7 y-1 
TL’ Ll ~ (P3/Pp) 
Recalling the definition of “nee 
: 1/2 lye (A-25) 
G = (1 + x) 
M “pp 
Mnese coefficients describe vpeecriormance ofan 
ejector as an energy or momentum producing device. Since 


an ejector is essentially a pump, it becomes necessary to 
develop certain coefficients whieh describe the pumping 
characteristics. 

The drivane pressurewcoet i semen: ae is defined as: 


P z= 
a ae Be (A22iay 


ia 


or, for an ideal gas, 


2 = ne) on val (Asam 
la ee ; _ 
Payt LE = F357?) G 


Ca 


D 








T 
aba 


The suction pressure coefficient is similarly defined. 


c= Pro 7 f3 (A-28) 
ae He 
OF 
, =(1 Mk 1 eee ( Ae 
ne |e i 
7 (P/P py? 


The geometric characteristics on an ejector may also 
be non-dimensionalized in the same manner. Consider the 


cross sectional area of the mixing section A , then 


5 
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ae 


C = a bi cane d (A-30) 
ee ro Tre My Me 


petting sup ene Tour dimensionammuniveequations yields, 





A, (L*) - aay 6? [e) [ Lag” (ESI) 
or 

eye 2 = -la -c + 2d (A-32) 

ee o = ate (A-33) 

T t= 2a +c = 2d Cho) 

Q: Con | (A-35) 


adding A-32 and A-34 gives 


eo = -2a 
a = -l 
aoe +h 
oe 
gd = -l/2 
chus 
-36 
Cag Ag wah fa, (A-36) 








a 
It was found that the actual controlling geometric 


parameters could be expressed by an area coefficient and the 


wae Oot As to the primary nozzle area Ay: nus | 
ee 
a (A-37) 


was used for the non-dimensional area parameter along with 


Ca3: 
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3. Detailed Development 
Consider the mixing Of Cwo cereemuric jets am a con 


Stant akheagewerricumonLecomalcius 






Primary 


Secondary a 
Nomenclature 
AREAS 
A, = Primary area 


> 
li 


5 secondary area 
Az = Mixing section area 


: ; : : 2 
introduce non-dimensional area ratio, K', thus, 


Hy 


Ao 


Se 
fe = K AS 


MASS FLOW RATES 


2 
(KX - 1) Ay 


ms = Primary mass flow rate 
m5 = Secondary flow rate 
m = Total flow rate 
ceca lini aeriniton Ol eine 2p 
a 
M5 = X mM 
DENSITIES 


p, = Primary density 


secondary density 


© 
ll 


EXivt density 


~ 
il 
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Lec 


Pa 
a= 
a 
VELOCITIES 
V, = Primary velocity at exit from nozzle 
V, = secondary velociy wae. teerirom WozZzZle 
V3 = Final velocity 
ree 
V 
r = 2 
7 
i 


From the definition of mass flow rate, 


m P5 Ao V 
xX = ae = === (A-38) 
eee 1 


Or finally. 
x = 6(K* -1)? (A-39) 


fet OM CG Olives yy. 


= my +m (A-40O) 


A, V, + p> Ag V 


ae 2 


Pp 
4 {a= (1 + x) 
ra he(¥2 
Chere eices 


meme tex) [Vy 
a K? (7 oa 


Now, consider the constant pressure mixing of two ideal 


gases ina constant area duct. 


os 





Cale + ms. Ceo 
"i, “pen "2 “p 92 eee 
or 
C.P, -xC_P, a 
a oo 
Py i PR 
S 
ee -l = tant 
ihOr abate = Sa 
RY 
therefore, 
P xP 
ae o = (1 pee 
i 2 
3 
uae Po — P, = Ps 
thus, 


| i to x7ce 


(1 + x) Se 


As a comparison, consider the mixing of two liquids 
at constant volume, 

1 +7 e 
oS 
P) Pp 

Or 
mee + x/6) 
—|= (A-43) 
P2 (1 + x) 

Now assume that the actual mixing process can be 
described in a similar manner such that the compressibima 
effects are taken into account by some change in volume of 
the primary fluid. By introducing a compressibility factor 


wu, the mixing process can bewadeeeribed by 
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_ tee a 
S 
The vyainewor WT will be umity for the case of incom- 
pressible mixing. For the compressible case, uw can be 
determined fromenergy considerations as will be shown later. 
Rewriting equation A-44 


Ee _ (1 + x) 
on Oe 7 0) 


Combining this result with equation A-41, 


. 





therefore, 


3] Git x/s) 
= z 


Ve i (A-45) 


al 


Now, apply the momentum equation for the case of a 


frictionless mixing section with uniform velocity distribu- 


tions. 
PL AL + P A, ~ Ps A. = m. Ve Fe my Vy - m, Us 
davidii ne Shiiepo y Py and Ay yields 
2 2 Dee 2 2 2 2 
Py + BD (KT-1) _ 2 a ("3 K°V3) - Vy - 6(K"-1)V,, 
eal Ell oa a 
(A-46) 


V 
after substituting for(®3| . ee 73] the following 
07 i. Vi 


expression is obtained: 


Ie P 2 2 
al D on _ P D (1 + x)(u + x/8) X San 
1 3 ee 2 | 2 BTK2-2) 7 : 

(A-47) 
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Nowe COlsmder 


Paso ee! 7 


1 
when 
Mai ei 
@ie 
Ps = Ps + AP. 


thus, looking at the left-hand side of equation A-47, 


2 
_ 2 
Pot AR, + Poh 2 eee 
p p 
Ps 1 1 
Or rearranges. 
Po 
2 AP 
Sasi nieon —--C 1- = 2 (A-48) 
0 
0,Vy iss 


Define a non-dimensional pressure rise across the 


Mixinee Secy beim. 





A P2- 
Pe re gt eg te (x) (tx) 
2a a §“(K2-1) sc Ke 
(Aes 


With this result, along with the definition of YW) fae 
possible to obtain an interaticneseiluvicn vo thesia. 
probDlemerer any given set Of Skievinieamemevers. 


Consider the following known parameters: 


es Ps = Pamb — fixed value 

O. Pry = Controllable value 

je Pmo = Controigete value 

4, To/Tpy = Temperature ratio can be controlled 
Dae Ke = Known fer any given ejecrcr 
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Now, for a specified value of each of these parameters, 
the following technique was employed to obtain a solution: 
a. Assume a value of M>- 


bee Hor M. < 1 jj. weaeutareo eli. 


1 1 

FE -1 

a = 1 + y-1 M Z ee 

Pl e al 
or 

iE -1 
na? =[7ér] fe] OF - 3 
i 

DUL 

Po 

ee 

eee ml J 
2 
therefore, for M, <i Po - Ps 
and 
—1 
2m. Pri| = y-1 ,, 2 
M =/{_c_ eae dips el eli -~50 
: ey 7 2 M2 soi: 

for M, >1, set M. = 1, since the primary nozzle 


oe i 


is convergent. 


c. Density Ratio 


en 22 eee hie 
Py (P17 Pp) Png 


One 
2 41/y-1 


y-l 
aa a — ae 
D 


Or 
ll 


(A-51) 


a Seo al M a bl Te 
2 
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qd. Pl@ee Rate Ratio 


mo _ P2 Ap Vo 


eee 
i eee 
ONG 
M 
eo (Koo) Li | Ts 
My pT 
1 
IDL 
oe (ef te) Te 
T,  (TY/T Hy) 
2 
= 2 Ie 1S i 


— +e 
+ aa “3M, 2 


there fore, 


See 
§(K2-1) Be | ze J 1+ em“) | 172 
Tl 


fie 
1+ O=n (A-52) 


e. From these inlet conditions, a value of 9% can 


be determined. 


P P 
a 3 ae: 2 
P> Vy 
(ah el (te 
ae a a a 
2 D 
aoe =e 
Py Pa 
coe = |{E3) ~ (PEN) / age? 
Pa es éyM, (A-53) 
look at 
P 
(RE) IE ots) 
a Poy} \ Py Py a 
1S +N, y/y-1 
eee © 2/0) a] - 2 +g 
pays p 
it Ca fF) VT 1 4 wees 


{0 


or finally, 


, 
d = Fa [2 i tla?) “= (fr 
iC Ee 2 TI 
y-1 y/y-1 
+ ae 2 
éyM 


14 ou, 1 (A-54) 
f. From the momentum definition of 6, and by 
assuming an maitial yaluc Of simeorsumiry. =a value of a can 


be determined. 


APS in x? 4 ee 


2 
Py VT ee sn §2 Ke 





ee 
ae 3 


(A~55) 


NOTE: All terms in this equation are known with 


the exception of AP} 








PoV 2 
al 
ec eae. = 0 
HOt see ale 
AP = Pl aan Po = 2528 Pm a P, 
and 
I 
x 2 
Po = ——-=—____-. y/y-1 
Cy + Yoly 7 
a 
nieise 1 Oren 
AP Pro Aes 
a 5 = E525 Pm ad 5 
p Y-ly, 2 
au {1 + ey | Poy 
or after some algebraic manipulation, 
-l y/yY-1 
ea iM, 23 | st “(52 ] 
Tg aie Pry 
y/y-1 
(A-56) 





ip ioe a is known. 

ge. The necessary condition for a solution is that 
e equal %,; thus, all boundary conditions are satisfied. 

ley 

ERF = 6, _ %; (A-57) 
such that for a given set of parameters and a particular 

Valle won sp, the SrOnee M5, has been chosen when ERF = 0. 

h. It is now possible to obtain a new value for uw. 


From energy considerations, 


Tomo = COne.r 


p 
Tp3 Bp (A-58) 
Tino 1+ x 
also, 
tee 
ie 3 ae 
p 
or 
T T a 
ak A ovR 
ne a y-l Tipe 


after some rearranging, 


(e)F2) [ped DT 


Ino Tipo 





from equation of state, 


so were [222 ae 


PT p T 
oe) 2° 3 2 Tia ae 
OY famaa Tile 
2 
0 V2 
oe o¥1° 4 o | { Tro 1 (A-60) 
— payor ear | |. oe 


tie 


Read A3Al1, GAMMA 






K= 1,7 
0 = 0.1(K) 


Do Loop 
P3PT1 = 1. 












Set 
PT2P1 = P3PT1 


U = 1.0 
ACH2 - 0.001 





M and PT2P1 
Iteration 


-.---4 Wega 1 


When two successive , 
values of M are _  , 
equal 


Calculate CPP, CPS, CM, CP, CS 


Print x, CPP, CPS, CM, CP, C 
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APPENDIX B 


Experimental Development 


DIcevne primary ObDjJCeCemve on mre wcxperamental program 
was to obtain a realistic comparison of actual ejector per- 
formance with the idealized theoretical predictions, the 
data obtained was reduced to the same non-dimensional 
performance coefficients as developed in the theoretical 
analysis. With the aid of a specifie computer program, 
EJECTOR, the raw data was converted to obtain these required 
eoerticients. 

1. Flow rate calculations 

As mentioned in Section 3, the primary and secondary 
flow rates were determined using squared edged orifice 
plates with flange pressure taps and an iron-constantan 
thermocouple in accordance with ASME standards. 


From Reference 4, 


m = =! Do° akyV ph (a 


m = mass flow rate 
Dye Oiiitee dilamever, 

a = thermal expansion factor 
Y, = compressibility factor 

K = discharge coefficient 

0 = f lade ncicy 


hw = pressure drop across the orifice 


fee 


Now, for an ideal gas, 


ae 


a|7 
i, 


with 
B it anches Ol snemeria,, 
LG alate OR 


hw in CM H50 


me 359.1 ( (Aaa 2 P_hw 
"* 3000 Vtes.sisyte.say 2 VA 








1 
; 9 P_ hw 
ioe(07205) De amanen 1 (B-2) 
e i 7 
1 
where 
Ty = static temperature Upectream Of 9the Ofna 


From Reference 5, for stainless steel orifice place. 
operating in a temperature range of 60° - 350° F. 


a = 1.0 + 0.0015 [Pare] 


007 (B-3) 
where 
a eo 460 
Cnheret ogee 
a = 1.0 + 1.5 (7 - 520) x 107? (B-4) 


Also, from Reference 5, the expansion factor 1oreaie 


Wit ey = eee 


Doe 
Y= Iso] (.041 = oon hw - 
, [041 35 (55 1 3 (B-5) 
Do 
for Nain Cleon HAO, P Sasi oe > = Dee? 
ilk 
= hw 
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The discharge coefficient K, was determined in order 


to correct for Reynolds Number effects. 


~ 6.316W 
NE Dp gl (B-7 ) 
2 
where 
W = 3600 m 
1 


4 = absolute fluid viscosity 
By determining a linear relationship between Zt and t 


from Reference 5, page 335, it was determined that 








ia aT 
42= 1004 =1. 0.24 
Gr 
2 =11.9 + 0.24 oy 6 | Ix Tose 
il 0 a Bae 
therefore: 
2 
Rp = | (6.316) (3600) x 10 (B-9) 


oT mn 
D,, 1.940 .24 [re - 5 .6| 
Now, for flange pressure taps the discharge coefficient 


K, can be determined 





| (B-10) 
R 
E 
where 
A =D, 1830 - 5000 | 22) 9000| 22 
} Dy dD 
1 
D 
- y200{_2 P 4 530 
Dy (B-11) 
VDT 
Ke = tabulated discharge coefficient for R, = ». 


Combining these results gives: 


Sal 


; 20 
= 00 Phw . 
See a (B-12) 


Now, for the primary air system 





C, = .07205 
Ke = .9756 
D, =e 067 in 
D, = 1.4465 in 

a S oa 
c t= C) Ke D, Oe, 

0 + LS aes 5.2| x 107° 

a = ; 2 100 . 
Y. = 1.0 - 0.0102 hw 

1 D 

h 
= 1 + 966.324 
Re 


Since ¢€ is dependent upon the actual fliow rate, the 
final solution is based upon an iteration solution, where 
originally ¢ is assumed to besumity> This Peeration ae 
performed with the CDC 1604 digital computer as is shown 
in the Fortran Program EJECTOR (Appendix B). The same 
procedure is followed in determining Phessecondary Vie 
rate. 

2. Average Exit Total Pressure 

In order to obtain the Stave vcencvien al Vice aap 
to the mixing section, a pitot-static Prandtl type preset 
probe was used to traverse the exit plane. The static pressure 
was found to be essentially constant as was the total 


temperature. The total pressure, however, varied along the 
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traverse. For this reason, it was necessary to determine 
a suitable equivalent averages vocal pressure. 

Consider a hypothetical constant area, frictionless 
duct appended to the actual mixing section. Assume that 
the velocity profile changes from the given non-uniform 
distribution to a uniform value. Since the Mach numbers 
aeeche exit plane were less than Ov 3s ecompressibility 


effects may be neglected. 





Actual Exit Hypothetical 
Plane Exit Plane 


For incompressible flow: 
S o 


and 
a = a 
Pm = Ee 1 / omer Vv (B-14) 


Pe Olas Corre inlet Vs 


fe VdAe= [ oVaGA 


+ 


A A 
Ore een ay cae 
V= | v dA (B-15) 
a A 


Now, applying the Momentum Equation, 


, 


jPod a-JP,aa=[oaan-[oveas (e=6) 
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fox Pa = const = 0 psig 
V a= Gomer 
oO = const 
Ce mmc OnSt 
thererore, 
ie us a - 
Py - P_ = -fov 1/2 » ¥ (Baies 
Due 
D_ Sm 72 fs 7 
T a q QO 
ye B 2g (Bailey 
Pp 
Now, substituting for V and V, 
q = [oe 2q dA i ae 
. K - 1/2 of [(2q dA : 
(=) (B-19) 
Cm polmee 0 as COMscanr 
oy ze 
q=2{ qaa- fal/? aa : 
A " A (B-20) 
however, since P amB = te ; 
H. = dA 1/2 
Hi 2 f hp SA - fo on (B-21) 


where Dn 


= total pressure head measure at the @€xXit pilamee 


This result was then converted to a finite difference 


format where: 


he 
A 


AA, 
A 


Hl 


total head measured at the z<th traverse 


annulus areas corresponding to the 4th traveree 
station 


there Cwewer 


2 
7 , 2 AA; 
Hq 22a as : 2 hye | (B-22) 
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After checking several calculations it was found that: 


2 
E nie mi = 2 ny ni (B-23) 


A A 


The approximation was then made that 
— = he? 
Hep cs h; = (B-24) 


This approximation gives an error of less than two 


percent for the worst case. 


Performance Coefficients 


Upon determining the individual flow rates and measur- 


ing the conditions at the exit plane, it was possible to 


calculate all performance coefficients. 


a. Pumping Power Coefficient 


pa ; 
Co 29Gn) far) 
PP =f ae 
ms H 
al eal? 
Cc. = (1+ x){_t3 8 
pp Tr Z 3 (B-25) 
i eas ‘ 
Pry of 
NOTE: Both Tp3 and a were measured directly 
a] 
along with Py. Png was determined from 


previous analysis (Section 2) while P. 


was ambient pressure. 





Oe, Suction Power Coefficient 
ne 
Cos = sae ie 
m4 ie 
y-1\(Tr2 = 
oe ae An Pr (B-26) 


Output Momentum Coefficient 


Cr A (my + mM ) V3 
m. Vi 
dl a 
eee + xX } cH3 
ae 
a 
- We lye 
Ou ee? oe C) ae) (B=24 


Driving Pressure Coefficient 


Pp. Ae 
_ Pry : Ps 
a | ye 
RImo} P Tl Pr 
Orweernally., Px 
ee: ie 
‘ps lfez|leay/ ah 
Ppa 1- (2 Nel 
aoe ly 
Pry 
SUCTION Pressure Coefificieny 
C = Poo a a 
2 Pmo H, 
Bere. 3) 
Pre a: Denies 
Rin |p Tl) To Y 
Poy 
Ores 
ep E 
=e ae om 


P Yeek 
a? 7894) 
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These coefficients were calculated from the experimental 


data with a rather basic data reduction program EJECTOR 


as Will be shown in the following section. 


4. Data Reduction Program -- EJECTOR 


The following is a comparison listing of Fortran 


Symbols to the actual Algebraic jymbols used throughout 


paas study: 


Algebraic 


Ke 


TARE 


Pamb 


Fortran 


A3Al 


TARE 


PAMB 


GAMMA 


H(T) 


ata 


2 


Ta) 
P(T) 
XP1 
XS1 


ACH2 


Paveinee 


2s 
Bees 


Remarks 


Ratio of the mixing section area 
to the primary nozzle area. 


Manometer atmospheric reference. 


Atmospheric pressure in inches of 
mercury. 


Ratio of specific heats. 


Recorded pressure readings at 
Stata eileen 


Primary temperature recorded 
upstream of the orifice plate. 


Secondary temperature recorded 
upstream of the orifice plate. 


Total temperature at Station I. 
Absolute pressure at Station I. 
Primary mass flow rate. 
Secondary mass flow rate. 


Secondary fluid mach number at 
the entrance to the mixing section. 


Ratio of the secondary-to-primary 
total pressures. 


Average exit total temperatures. 


Average exit total pressure. 
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Algebraic 


Xx 


Fortran 


X 


Cr 


CPS 


CP 


CS 


CM 


Remarks 


Ratio of secondary-to-primnary mass 
flow rates. 


Pumping power coefficienc. 
suction power coerricient. 
Driving pressure coefficient. 
SUCTION pwesstece cocrii¢tenie 


Momentum coefficient. 


The program, EJECTOR, was set Up to read in any numa 


of sets of data. 


For each setj;~a@ new value of A3Al1, PANES 


TARE, and GAMMA must be read. The program consists of two 


overall do loops combining the two flow rate iterations =a 


the calculation of all required coefficients. The lorie 


flow chart and Fortran program are shown on the following 


pages. 
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READ VALUE OF NSETS 


READ TARE, A34 PAMB, GAMMA 







Print Heading 
for 


Print Out 


DO LOOP I = 1, NSETS 


READ H(L) ) 


x f" A] D 


DO LOOP to Convert 
Gage Pressures to 
ABSOLUTE 
P(I) = TARE + PAMB-H(T) 






Flow Rate Iterations 


XSO = 

XS1 = XSO 

ZET1 = 

XS2 = XSO*ZET1 


Compare XSI & XS2 


M2 Iteration 
Calculate F 
M21 = F 


M22 = F(1 + .2 M21 ) 
Compare M21 & M22 





Calculate CPP, CPS, CM, CP 


Print Out 
GEP, CPay Ch CrP 
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